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ABSTRACT 


Spatial  distributions  of  the  ground  clutter  backscatter  coefficient  o°  for 
various  types  of  terrain  have  been  found  to  fit  a  Weibull  distribution  quite  well. 
This  report  discusses  the  Weibull  probability  distribution  function  and  its  char¬ 
acteristics,  and  presents  some  clutter  measurements  that  demonstrate  the 
Weibull  distribution.  The  data  presented  include  measurements  taken  at  L-, 

S-,  and  X-band  frequencies,  at  several  depression  angles,  and  for  various 
resolution  cell  sizes.  A  method  for  generation  of  random  variables  that  fit  a 
selected  Weibull  statistical  population  for  use  in  radar  simulations  is  also  given. 
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a 

b,  o 


P 


r 


i 


a 

m 


o 


cr 


slope  parameter  (equal  to  l/b) 

Weibull  distribution  parameters 

probability  that  cr  0  ^  a0 
c 

,th 

the  1  random  number  of  a  uniform  distribution 

clutter  backscatter  coefficient 

mean  value  of  cr° 

cr°  for  any  particular  clutter  cell 

median  value  of  a0 
th 

the  i  value  of  o°  from  a  Weibull  population 
variance  of  the  Weibull  function 


IV 


1.  Introduction 

Procedures  for  estimating  the  performance  of  a  radar  in  a  clutter 
environment  vary  widely.  Probably  the  most  common  approach  is  to  treat 
ground  clutter  as  though  it  were  homogeneous  and  characterized  by  a  constant 
hackscattcr  coefficient,  tr°.  Depending  on  the  particular  value  chosen  for  a0, 
the  detection  range  can  vary  widely  since  the  stgnal-to-clutter  ratio  is  generally 
considered  to  vary  inversely  with  range  to  the  first  power. 

Other  theoretical  models  utilize  a  statistical  clutter  mode!  L  describe 
the  spatial  distribution  of  land  clutter.  The  most  common  probability  distribu¬ 
tions  used  to  model  the  spatial  distribution  of  a0  are  the  log-normal  and  expo¬ 
nential  (or  liayleigh  amplitude)  probability  density  functions.  At  low  depres¬ 
sion  angles,  the  land  clutter  density  generally  varies  over  a  wider  dynamic 
range  than  would  be  predicted  by  the  Rayleigh  amplitude  distribution;  and  for 
this  reason,  the  log-normal  distribution  is  more  widely  used  for  low  depression 
angles  to  model  land  clutter  spatially. 

Airborne  land  clutter  measurements  made  by  the  Naval  Research 
Laboratory  [1]  were  found  to  fit  approximately  the  Rayleigh  amplitude  distribu¬ 
tion  in  a  number  of  cases.  These  measurements,  however,  were  made  at 
depression  angles  on  the  order  of  5  degrees  or  more,  which  may  have  accounted 
for  this  type  of  distribution. 

Another  probability  function,  the  Weibull  distribution  function,  has  been 
observed  to  fit  some  types  of  land  clutter  spatial  distributions  of  cr°  quite  well. 
The  Rayleigh  amplitude  probability  distribution  function  is  a  special  case  of  the 
Weibull  function,  and  preliminary  results  appear  to  indicate  that  the  clutter  dis¬ 
tributions  approach  the  Rayleigh  case  for  larger  depression  angles. 

This  report  discusses  the  Weibull  probability  distribution  function  and  its 
characteristics,  and  presents  some  clutter  measurements  that  demonstrate  the 
Weibull  distribution. 


2.  The  Weibull  Distribution  Applied  to 

The  Weibull  probability  density  function  is  a  single  variate  function 
having  two  parameters,  a  and  b,  and  is  given  by 

p(<r°)  =  — ^ - exp^-(7ob/o^  ,  <1 


\ 


i 


where  a0  is  the  variate  in  terms  of  the  clutter  backseatler  coefficient.  I'or  the 
case  where  b  =  1,  Kquation  (1)  becomes  the  exponential  (or  Rayleigh  amplitude) 
density  function 


.0  1  ,  0  ■'  0. 
p(n  j  = - exp  t-o"/CT  )  , 


(?) 


where  o  =  cr',  the  average  value  of  cr°. 

Integration  of  Equation  (1)  yields  the  probability  that  the  actual  clutter 
cross  section  per  unit  area  will  not  exceed  some  value  rr°.  resulting  in 


P  =  1  -  exp  ^-a°b/o^ 


<:i) 


The  mean  of  the  Wcihull  distribution  function  may  be  determined  from 
the  integral 


r°  =  f  tr°p((T0)  dn0  . 


H) 


A  eiur  ge  of  vi  riyble  x  =  oob/o  (and  defining  a  change  of  parameter  a  =  l/b) 
yields  the  it  tegra' 


~o  a  r  a  -x 
a  =  a  J  x  e  dx  , 


(5) 


which  is  the  gamma  function  integral.  The  resulting  mean  value  of  n0  is  then 

<r°  =  fyar  ( 1  +  a)  .  (b) 

Equation  (d)  v.an  be  solved  for  <7°  to  obtain 


a  t 


n(T b) 


'  l/b 

/  i  \ 

- 

a  tn( - 

\  1  "  P/ 

(7) 


The  m?di:-r  of  the  distribution  is  then  determined  by  setting  P  =  1/2  in 
Equation  which  reduces  to 


a  =  fa  fn  2] 
m 


(8) 


The  parameter  c*  can  nov.  be  expressed  in  terms  of  the  median  clutter  back- 

scatter  coefficient  a  0 
m 


a- 


fn  2 


(0) 


Substitution  of  Equation  (9)  in  Equation  (3)  results  in  the  Weibul!  distribution 
function  expressed  as 

P ^crc°  -  =  1  -  exp^-fn  2  (cr%m0^  U]  .  (10) 


Analytical  techniques  are  available  for  determining  how  well  a  measured 
distribution  of  the  clutter  backscatter  coefficient  approaches  a  Weibull  distribu¬ 
tion  and  for  estimating  the  parameters  "b"  and  (or  a  if  desired) .  A  graph¬ 
ical  approach,  however,  is  useful  and  more  easily  applied  to  experimental  data, 
and  this  method  rather  than  the  analytical  approach  will  be  described  here. 


Substituting  the  results  of  Equation  (9)  into  Equation  (7)  ,  one  has 

(in  2)  v  ' 


(ID 


Taking  the  logarithm  of  o°  and  multiplying  by  ten  results  in 

10  log  cr°  =  10  log  -  10  a  log  ( ?  n  2)  +  10  a  log  t~  ^ — * 


n(rT>)]  • 


(12) 


which  is 


cru(dB)  =  <7ni°(dB)  +  1.6  a  +  10  a  log 


(13) 


Now  it  car.  be  seen  that  if  the  Weibull  probability  distribution  is  plotted  on  graph 


paper  having  the  probability  scale  proportional  to  log 


and  a°  in 


decibels  on  a  linear  scale,  the  result  is  a  straight  line  hating  a  slope  determined 
by  the  parameter  "a.  "  The  value  of  the  median  fixes  the  relative  position  of  the 
distribution  line.  The  graph  paper  having  these  special  scales  described  above 
is  sometimes  referred  to  as  Rayleigh  paper,  since  the  Rayleigh  amplitude  dis¬ 
tribution  is  a  special  case  of  the  Weibull  type  distribution  for  which  a  =  1 . 
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Measured  distributions  of  a0  can  be  plotted  on  Wcibutl  paper  to  determines 
graphically  if  the  distribution  is  of  the  Weibull  type,  and  if  so,  to  determine  the 
value  of  "a"  (or  b) .  If  the  points  lie  in  a  reasonably  straight  line,  the  Weibull 
distribution  can  be  considered  to  be  a  reasonable  approximation  to  the  actual  dis¬ 
tribution,  especially  if  the  points  fit  well  within  the  10-  to  90-percent  range. 
Values  of  a°  corresponding  to  the  lower  probabilities  are  generally  well  below  a 
system's  noise  level  and  present  virtually  no  problem  to  the  radar,  while  strong 
values  of  cr°  corresponding  to  the  higher  probabilities  arc  generally  too  large  to 
be  suppressed  effectively  by  moving  target  indicator  (MTl)  systems,  anyway. 

The  capability  of  a  radar  to  operate  satisfactorily  in  a  land  clutter  environment 
is  strongly  determined  by  the  manner  in  which  the  clutter  is  distributed  spatially 
over  the  range  of  cr°  which  might  be  described  as  moderate.  The  use  of  any  one 
particular  vaiue  of  o0  by  itself  tends  to  be  misleading  in  describing  a  radar's 
performance  from  a  system  standpoint. 


3.  Generation  of  Weibull  Distributions  in  Simulotions 

The  Weibull  distribution  function  is  readily  adaptable  for  use  in  digi¬ 
tal  computer  simulations  that  utilize  a  Monte  Carlo  approach  to  simulate  the  per¬ 
formance  of  radar  systems  in  a  ground  clutter  environment.  As  a  target  moves 
over  the  ground,  it  competes  with  different  clutter  cells  or  with  different  com¬ 
binations  of  clutter .  If  it  can  be  assumed  that  there  is  little  or  no  correlation 
between  successive  looks  of  the  radar  at  a  particular  target,  Equation  (11)  can 
be  used  with  a  random  number  generator  to  produce  random  variates  <r  °'s  whose 

density  function  is  Weibull.  The  above  assumption  is  generally  valid  for  search 
radars  with  good  spatial  resolution  and  frame  times  on  the  order  of  several 
seconds  against  tactical  aircraft. 

Equation  (11)  is  simply  an  expression  which  will  transform  or  map  a  set 
of  uniformly  distributed  random  numbers  which  cover  the  unit  interval  between 
0  and  1  into  the  domain  of  a0.  To  generate  random  cr0,s  from  a  statistical  popula¬ 
tion  whose  distribution  is  Weibull,  uniformly  distributed  random  numbers 
between  0  and  1  are  first  generated.  Each  random  number  r.  is  then  transformed 
into  a  corresponding  <r.°  by 


a 


o 

i 


(fn  2)a 


(14) 


where  the  slope  parameter  "a"  and  the  median  o^0  have  been  previously  defined 
for  the  particular  type  of  clutter  to  be  simulated. 
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The  Weibull  distribution  function  has  been  used  to  simulate  spatial  ground 

clutter  distributions  in  the  lb  S.  Army  Missile  Command  Air  Defense  Simulation 

(MADS)  digital  computer  model.  Tills  program  utilizes  digital  terrain  data  at 

each  500-meter  interval  over  the  area  simulated  for  use  in  detei  mining  masking 

effects.  A  code  number  is  used  to  identify  the  type  of  ground  cluttei  (forest, 

grassland,  etc. )  at  each  500-meter  grid  coordinate.  The  code  number  selects 

the  desired  Weibull  distribution  having  slope  parameter  "a"  and  median  a  0  for 

m 

a  particular  clutter  type,  and  a  value  of  a0  is  drawn  from  the  selected  Weibull 
population. 


4.  Spatial  Backscatter  Distributions  of  Measured  Lana  Clutter 

Several  sources  of  measured  distributions  of  clutter  backscatter 
coefficient  have  been  obtained  which  exhibit  the  Weibull  distribution.  The  data 
points  have  been  plotted  on  Weibull  paper  as  described  in  Part  2  such  that  a  dis¬ 
tribution  which  is  Weibull  will  result  in  a  straight  line.  These  sources  represent 
measurements  from  both  high  and  low  resolution  radars,  relatively  long  and 
short  clutter  cell  averaging  times,  and  for  various  frequencies.  The  greatest 
deviations  from  the  straight  line  appear  at  the  very  small  values  of  a0  which 
might  have  been  a  result  of  the  receiver  noise  level. 

Figure  1  shows  a  distribution  of  Rocky  Mountain  [2-4]  ground  clutter 
taken  at  S-band  with  a  resolution  cell  defined  by  a  2.  0-microsecond  pulse  and  a 
1. 5-degree  beaniwidth.  The  data  were  taken  for  ranges  less  than  20  miles  from 
an  unmasked  location.  No  information  is  giver  concerning  the  depression  angles 
or  vegetation.  Nevertheless,  the  distribution  fits  a  straight  line  quite  well  on 
Weibull  paper  over  a  wide  portion  of  the  probability  scale.  The  circled  points 
indicate  values  taken  from  the  distribution  curve  previously  given  [2] . 

The  distribution  of  the  clutter  backscatter  coefficient  for  data  taken  near 
Huntsville,  Alabama,  |5]  is  shown  in  Figure  2.  These  measurements  were 
made  at  L-band  out  to  a  maximum  range  of  approximately  40  kilometers  with  a 
resolution  cell  defined  by  a  3.0-microsecond  pulse  length  and  a  1. 7-degree 
beaniwidth.  The  terrain  included  grassland,  low  rolling  wooded  hills,  wooded 
mountains,  and  some  man-made  structures.  The  data,  taken  with  a  nonscanning 
antenna,  and  an  averaging  time  of  3  seconds  at  each  position,  represent  the 
strongest  241  clutter  cells  out  of  several  thousand  cells  examined  and  therefore 
exhibit  a  higher  median  value  when  compared  to  other  measurements.  The  radar 
site  was  located  on  one  of  the  hills  such  that  the  depression  angle  from  the 
antenna,  as  detet  mined  Irom  maps  of  the  area,  was  generally  less  than 
0.5  degree.  Even  though  the  smaller  clutter  cells  were  neglected,  the  Weibull 
distribution  still  represents  a  good  fit  to  the  distribution  of  the  measured  clutter 
data. 


Data  taken  at  X-band  by  Linell  [6]  at  the  Research  Institute  of  National 
Defense  in  Sweden  is  presented  in  Figures  .‘i  through  5.  The  resolution  ce'l  is 
defined  by  a  0. 17-microsecond  pulse  length  and  a  1. 4-degree  beamwidth.  The 
data  were  taken  with  a  rotating  antenna  atop  a  water  tower.  Figure  3  shows  two 
spatial  distributions  of  the  clutter  backscatter  coefficient  for  a  Swedish  forest 
area  at  different  times  of  year.  These  particular  data  are  for  a  depression 
angle  of  0.7  degree.  Seasonal  deviations  appeared  only  during  the  month  of 
November,  resulting  in  an  apparent  shift  in  the  distribution  of  about  6  decibels. 
The  slope  changed  only  slightly. 

Figures  4  and  5-show  spatial  distributions  for  cultivated  land  taken  in 
the  months  of  April  and  May  at  varying  depression  angles  of  1.25,  2.5,  and 
5.0  degrees.  As  the  depression  angle  ^and  consequently  the  median  uj  ^ 

increases,  the  slope  of  the  distribution  approaches  that  of  the  exponential  (or 
Rayleigh  amplitude)  distribution.  This  trend  was  also  evident  in  measurements 
of  the  same  terrain  at  other  seasons  of  the  year.  Naval  Research  Laboratory 
airborne  measurements  [1]  at  several  frequencies  where  depression  angles 
were  on  the  order  of  5  degrees  or  greater  also  approached  the  Rayleigh  ampli¬ 
tude  distribution.  This  effect  is  probably  attributed  to  increased  shadowing 
effects  at  smaller  depression  angles.  Hopefully,  satisfactory  relationships 
between  the  slope  parameter  ''a,  "  the  median  a  °,  the  frequency,  the  polariza¬ 
tion,  and  the  depression  angle  (or  grazing  angle)  for  various  terrain  types  can 
be  verified  for  use  with  Equation  (14)  in  a  spatial  clutter  model  for  air  defense 
simulations. 

The  parameters  for  the  distributions  of  this  section  are  given  in  Table  I. 

The  median  a  l;  is  obtained  by  simply  reading  it  from  the  distribution  at 
m  ' 

P  =  50  percent.  The  slope  parameter  "a”  is  obtained  by  .he  method  described 
in  Appendix  B. 
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PROBABILITY  THAT  CLUTTER  DENSITY 


PROBABILITY  TNA 


PROBABILITY  THAT  CLUTTER  DENSITY  DOES  NOT  EXCEED  o°  (p*rc*nt) 


FIGURE  3.  GIwTTND  CLUTTER  SPATIAL  DISTRIBUTIONS 
FOR  FOREST  AT  DIFFERENT  TIMES  OF  YEAR 
AT  X-BAND 
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PROBABILITY  THAT  CLUTTER  DENSITY  DOES  HOT  EXCEED  «° 


FIGURE  4.  GROUND  CLUTTER  SPATIAL  DISTRIBUTIONS 
FOR  CULTIVATED  LAND  IN  APRIL  AT  SEVERAL 
DEPRESSION  ANGLES 
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FIGURE  5.  GROUND  CLUTTER  SPATIAL  DISTRIBUTIONS 
FOR  CULTIVATED  LAND  IN  MAY  AT  SEVERAL 
DEPRESSION  ANGLES 
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Appendix  A 

DETERMINATION  OF  THE  VARIANCE  OF  THE 
WEIBULL  FUNCTION 

In  determining  the  variance,  the  mean  square  value  of  u°  is  first 
calculated: 

~  2  "  2  Vy.0^  1  /  b/  \ 

<j°  =  f  cr°  p(c°)dc°  =  fa0 -  exp  (-a0/ a)  da0  .  (A-l) 

j  j  a 

0  0 

b/ 

A  change  of  variable  x  =  cry  a  (and  defining  a  change  of  parameter  a  =  l/b) 
yields  the  integral 

0  f  2a  2a  -x 
cj  -  J  a  x  e  dx 

0 

'-a 

<r°  =  o: *"  T(1  +  2a) 

From  equation  (6) ,  the  mean  is  shown  to  be 

<r°  -  oar(l  +  a)  ;  (A-4) 

the  variance  is  found  from 

a2  =  <7°'  -  (^)2  .  (A-5) 

The  variance  is  then 

a2  =  a2a[r<l  +  2a)  -  r2(i  +  a)  ]  .  (A-C) 


(A-2) 

(A -3) 
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Appendix  8 

DETERMINATION  OF  SLOPE  PARAMETER  "a" 


The  slope  parameter  "a"  of  a  Weibull  distribution  may  be  determined 
directly  from  che  plot  of  the  distribution  on  Weibull  paper.  By  recalling  that  the 
probability  scale  is  proportional  to  log  [f  n  (l/l  -  P)  ] ,  a  variable  is  introduced 
x  =  log  [fn(l/l  -  P)  ] ,  such  that  the  probability  scale  is  transformed  into  linear 
units  of  x.  Substituting  x  into  Equation  (13)  gives 

cr°(dB)  =  CTm°(dB)  +  1-6  a  +  10  ax  .  (B-l) 

The  desired  slope  is  the  derivative  — 7  : 

dx 


du°(dB)  _  _  A<r°(dB) 

dx  a  Ax 

The  slope  parameter  "a"  is  thenTound  from 

=  Ag°<dB) 
a  10  Ax 


(B-2) 


(B-3) 


It  is  now  necessary  to  determine  what  interval  of  probability  corresponds 
to  a  Ax  of  unity.  The  lowest  probability  shown  on  the  vYeibull  paper  used  in  this 
report  is  P  =  0.10.  Thus,  for  P  =  0.10, 


xr-o.i  ■  ,og 


['■tt] 


=  log  0. 1053 


xp=o.i  =  -°-9775  • 


Then,  for  unity  Ax, 


A  x  =  1  = 


P=0. 1 


and 


x' 


1  +  X 


P=  0. 1 


0.0225  . 
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Solving  the-  expression 


log  (fn  (1/1  -  P')  ]  =  0.0225 

for  P',  one  obtains  the  value  P'  =  0.051.  Thus,  a  Ax  of  unity  corresponds  to 
the  interval  on  the  probability  scale  between  10  and  f>o.  1  percent.  From  the 
above  result,  it  can  be  seen  that  the  slope  parameter  "a"  is  found  by  the 
expression 


a 


g(dB)P=0.G51  -*(dDWl 

10 


(B— 1) 


Graphically,  this  is  easily  determined  by  translating  the  slope  of  the  distribution 

over  to  the  probability  scale  such  that  it  passes  through  P  =  0.G51  on  the  right- 

hand  probability  axis.  The  interval  of  cr°(dB)  between  the  right-hand  probability 

axis  and  the  line's  intersection  on  the  P  =  0. 10  axis  is  /  <x°,v  ,  -  a0,,  „  \ 

\  P  =  0.  Gol  P=0. 1 ) 

above.  This  graphical  technique  is  illustrated  for  several  slopes  in  Figure  B-l. 
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PROBABILITY  THAT  CLUTTER  DENSITY  DOES  NOT  EXCEED  «°  (p.rcnt) 


FIGURE  B-l.  EXAMPLE  SLOPES  FOR  WTIBULL  DISTRIBUTIONS 
HAVING  DIFFERENT  SLOPE  PARAMETERS  "a" 
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